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Sperm cooperation has evolved in a variety of taxa and is often considered 
a response to sperm competition, yet the benefit of this form of collective 
movement remains unclear. Here, we use fine-scale imaging and a mini- 
mal mathematical model to study sperm aggregation in the rodent genus 
Peromyscus. We demonstrate that as the number of sperm cells in an aggregate 
increase, the group moves with more persistent linearity but without increas- 
ing speed. This benefit, however, is offset in larger aggregates as the geometry 
of the group forces sperm to swim against one another. The result is a non- 
monotonic relationship between aggregate size and average velocity with 
both a theoretically predicted and empirically observed optimum of six to 
seven sperm per aggregate. To understand the role of sexual selection in driv- 
ing these sperm group dynamics, we compared two sister-species with 
divergent mating systems. We find that sperm of Peromyscus maniculatus 
(highly promiscuous), which have evolved under intense competition, form 
optimal-sized aggregates more often than sperm of Peromyscus polionotus 
(strictly monogamous), which lack competition. Our combined mathematical 
and experimental study of coordinated sperm movement reveals the impor- 
tance of geometry, motion and group size on sperm velocity and suggests 
how these physical variables interact with evolutionary selective pressures 
to regulate cooperation in competitive environments. 



1. Introduction 

The factors that contribute to reproductive success are numerous and complex, 
yet across vertebrates, relative sperm motility is often the best predictor of male 
fertility [1-7]. When competition among males intensifies, adaptations that 
improve sperm swimming performance are therefore expected to be strongly 
favoured [8,9]. Indeed, comparisons between related taxa reveal that sperm 
of polyandrous species, in which females mate with multiple partners during 
a reproductive cycle, swim faster than sperm from closely related monogamous 
species [10,11]. Among the many strategies that improve sperm swimming 
performance, perhaps the most intriguing mechanism involves cooperation or 
association with other motile cells [12]. Even without direct attachment, 
sperm of some species interact with one another via flow fields that result 
from hydrodynamic interactions [13]. These associations, however, are magni- 
fied when multicellular groups form by conjugation, ranging in size from 
sperm pairs to large aggregates containing hundreds of sperm (reviewed in 
[14,15]). Sperm aggregation is often assumed to improve motility, yet compara- 
tive studies have shown inconsistent results (reviewed in [14,15]), and the 
underlying mechanics of the associations remain largely unknown. 

Like most muroid rodents, sperm from mice in the genus Peromyscus typically 
possess an apical hook on the head (figure la-c) [16] that is thought to facilitate 
the formation [17] and/ or stabilization [18] of sperm groups (but see [19]). Aggre- 
gations of Peromyscus sperm cells are formed by secondary conjugation [12]: 
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Figure 1. Scanning electron micrographs of [a) whole Peromyscus sperm cells, 
and head morphology of a single (b) P. maniculatus and (c) R polionotus sperm, 
(cf) Head orientation of sperm in a typical aggregate with hooks facing inward, 
and aggregates consisting of (e) two, ( f ) seven and (g) thirteen P maniculatus 
cells. 



sperm are ejaculated as solitary cells, but quickly begin to 
form multicellular aggregates by adhering to one another at 
or near the hook (figure Id) [20]. Overall, these motile sperm 
groups have a larger average velocity (straightline velocity, 
VSL; figure 2) than single cells; however, the largest groups, 
those over twenty cells, are often not motile at all [20]. Under- 
standing how sperm aggregates achieve greater average 
velocity than single cells, whether by increasing their speed 
(curvilinear velocity, VCL; figure 2) or travelling in a straighter 
trajectory (linearity), and how group size can hinder motility, is 
key to understanding how post-copulatory male -male compe- 
tition may be acting on sperm behaviour to drive and constrain 
group formation. 

In the genus Peromyscus, sperm competition is predicted to be 
greatest in P. maniculatus, because both sexes mate with multiple 
partners, often in overlapping series just minutes apart [21], and 
females frequently carry multiple-paternity litters in the wild 
[22]. By contrast, its sister species, Peromyscus polionotus, is strictly 
monogamous on the basis of both behavioural [23] and genetic 
data [24]. The sperm of both species form aggregations with simi- 
lar geometry and cell orientation, probably owing to analogous 
morphology of their sperm heads [25], yet the competitive 
environments experienced by P. maniculatus and P. polionotus 
sperm represent divergent selective regimes, which is believed 
to shape how cooperative sperm groups assemble [20]. Here, 
we use a minimal mathematical model to predict how sperm 
can improve their average velocity by forming aggregations 
and then use fine-scale imaging to test these predictions and 
gain a deeper understanding of how sexual selection has acted 
on this unique form of cooperation in Peromyscus sperm. 

2. Material and methods 
(a) Mathematical model 

A simple mechanistic picture of how the average velocity of 
sperm is a non-monotonic function of aggregate size is suggested 
by the geometry of the aggregates shown in figure le-g. As 
sperm cells form small oriented clusters, their motive force and 




Figure 2. Schematic of the average velocity (VSL) and speed (VCL). VSL is 
calculated by dividing the distance between the initial and final position 
in a sperm trajectory (dashed line) by the time Af employed to move; 
VCL is found by dividing the length of the actual curvilinear trajectory 
(solid line) by Af. (Online version in colour.) 

cluster geometry can increase owing to the collective beating 
of their flagella that leads to a greater dynamical persistence. 
However, in large clusters, the geometry of the aggregate 
approaches that of an isotropic cluster so that their collective abil- 
ity to move is severely hindered. A minimal model described 
below allows us to quantify the advantage of cooperation in a 
competitive environment using observable physical variables. 

Our approach follows a set of models originally developed for 
flocking behaviour of organisms [26,27], which have been used 
successfully to describe collective motion in a variety of natural 
and artificial systems, including fish and birds [28], insects [29], 
bacterial colonies [30] and robots [31] (for details, see appendix). 
In this spirit, we treat sperm as individual self-propelled parti- 
cles [32] that can interact with each other geometrically and 
mechanically, consistent with the biology of Peromyscus sperm 
aggregation [20]. We restrict our attention to the d3rnamics of the 
aggregates once they form, not attempting to address the process 
of hydrod3rnamic self-organization itself. Our method relies on 
three basic assumptions: (i) although the flagellum is responsible 
for propulsion, it does not contribute to mechanical interactions 
between sperm; (ii) the main physical mechanism associated 
with aggregate formation is due to adhesion between sperm 
heads, consistent with our understanding of sperm morphology 
[12,20]; and (iii) hydrodynamic interactions between sperm 
in the aggregate are negligible. Thus, although hydrod3rnamic 
interactions among neighbouring sperm are important in creat- 
ing self-organized patterns of swimming [33-35], in our minimal 
model that focuses on the dynamics of the aggregate, these 
interactions do not play a critical role. 

With the aim of characterizing the empirical system using a 
small number of experimentally measurable parameters, we con- 
sider exclusively those features of sperm mechanics that are 
essential for the formation of motile aggregates. Thus, we note 
that individual sperm occupy space, are able to move and can 
link to other sperm. Single sperm cells are then represented as 
two-dimensional tailless elliptical particles that self-propel at con- 
stant velocity Vq in a plane in the direction of their major axis n 
while being subjected to random planar rotations. Each particle 
is assumed to have a given number of 'keys' and Tocks', represent- 
ing the adhesion complexes on the sperm head. When the key of a 
particle is within a certain distance r^ from the lock of another par- 
ticles, a link, represented by a linear spring of stiffness /Ca, is formed 
(figure ?>a). If the key -lock distance eventually becomes larger 
than fa, the link breaks and the two sperm unbind (i.e. an indivi- 
dual adhesion complex can withstand forces up to a stall force 
^^a = I^Ja}- Finally, the particles are themselves assumed to be 
hard and unable to overlap, so that when in contact they pack as 
dictated by their geometry. 




(d) aspect ratio (a/b) (e) scaled adhesion range (rja) (/) scaled propulsion velocity iy^Tja) 
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Figure 3. [a) Schematic of the adhesive interactions modelled in equation (2.2). Sperm heads are treated as self-propelled elliptical particles whose major and minor 
semiaxes have length a and b, respectively. Each partide is equipped with a given numbers of keys and locks, representing the adhesion complexes where the sperm can 
bind. When the key of a partide is within a certain distance from the lock of another partide, a link represented by a linear spring is formed. The geometry of the 
aggregates affects their motility, so asymmetric aggregates [b] move fast and maintain a straight trajectory, whereas star-shaped aggregates (c) move slowly, because the 
velocities of the individual cells in the aggregate cancel each other. Average velodty [d-f) and speed [g-i) versus aggregate size obtained from a numerical integration of 
equations (2.1)-(2.5) for various aspect ratios alb (d,g), scaled adhesion range rja (e,h), expressed in units of the partide major semiaxis length a, and the scaled 
propulsion velodty VqTc/q (with = MD the timescale assodated with the rotational noise) (fj). (Online version in colour.) 



The above-described behaviour leads to equations of motion 
for the position of the fth sperm given by r/(f ) and its orientation 
di{t) relative to the x-axis of the laboratory frame given by 

^ = Voni + ^jL-'Fi and ^ = + (2.1) 

where the fth sperm has its major axis along W/ = (cos di, sin di), 
Fi is the total force acting on the fth particle resulting from the 
short-range steric interactions w^ith the neighbours and adhesion: 

Ni Li 

and Mi is total torque acting on the fth particle: 

M,= J2(dvxfrj)-l (2.3) 



v\^here N/ is the number of neighbours of the fth particle, L/ is the 
number of adhesive links, ks is the elastic constant associated 
with steric interactions, is the adhesive spring elastic constant, 
with ks ^ ka. Furthermore, 5/y and £ij represent the length of the 
springs associated with the steric and adhesive interactions, with 
Njj and unit vectors in the direction of the springs, /ul and y are 
translational and rotational drag coefficients, fij is any of the 
force between the /th and ;th cell appearing in equation (2.2) 
and dij the associated lever arm, z is the normal to the 
two-dimensional plane of motion and is a zero-mean 
delta-correlated Gaussian random variable: 

mm = 2DS{t-t'), (2.4) 

where D is a rotational diffusion coefficient. Here, we assume 
that the motion of the sperm is inertialess, consistent with the 
low Reynolds number regime they operate in, and further have 



ignored the effect of randomness in the translational degrees of 
freedom for the sake of simplicity. 

Our minimal mechanistic model of interacting sperm cap- 
tures the geometry of the individual sperm, their autonomous 
movement and finally their ability to interact with each other 
adhesively without overlap. While there are many possible var- 
iants of these models, the critical parameters in all of them will 
be qualitatively similar: the aspect ratio of the sperm head, the 
scaled ratio of the rotational Brownian motion to the interaction 
torque between cells, the scaled ratio of the adhesive bond 
strength to random fluctuations and the relative orientation of 
the adhesive bonds. These parameters together characterize the 
dynamics and persistence of movement in aggregates. 



(b) Sperm imaging and analysis 

Captive stocks of wild-derived Peromyscus maniculatus bairdii and 
Peromyscus polionotus subgriseus were obtained originally from 
the Peromyscus Genetic Stock Center and have been maintained 
at the Harvard University in accordance with guidelines estab- 
lished by Harvard's Institutional Animal Care and Use 
Committee. We used adult (age more than 90 days) sexually 
mature P. polionotus (n = 9) and P. maniculatus (n = 9) males for 
cross-species comparisons. 

After sacrifice, we immediately removed the left caudal epi- 
didymis of each male, made a single small incision in the 
tissue, submersed it in 1 ml of warmed Diggers -Whitten- 
Whittingham media [36], and incubated the tissue for 10 min 
at 37°C to release motile sperm. After the 10 min incubation, 
we removed the epididymal tissue, gently swirled the media 
and incubated for another 5 min. We collected 20 julI of media 
containing live sperm just below the surface of the aliquot, 
to reduce the number of dead cells, which sink to the bottom. 
We placed the aliquot on a plastic microscope slide and covered 
the sample with a plastic coverslip (plastic reduces adhesion of 
sperm to the slide compared with glass products), and recorded 
three 5 s videos (30 frames per second) of live sperm at 100 x 
magnification under phase contrast conditions on an upright 
microscope (AxioImager.Al, Zeiss, Jena, Germany). 

To examine the dynamic performance of sperm aggregates, 
we quantified the speed and velocity of both single cells and 
aggregated groups. The speed, also referred to as VCL, charac- 
terizes the rate of change of the two-dimensional projection of 
an aggregate's trajectory over time (figure 2). The average vel- 
ocity, or VSL, is defined as the rate of change of the projected 
distance along the vector connecting the initial and final point 
in the trajectory (figure 2). We acquired VSL and VCL data 
from video using the computer-assisted sperm analyser plugin 
for NIH ImageJ [37], which tracks motile sperm cells or groups 
to calculate VSL and VCL. We then estimated average 
linearity (VSL /VCL) for each track. Specifically, for each video 
recorder, we first used the 'find edges' and 'threshold' func- 
tions to isolate sperm images from the background and 
imposed a filter to discard tracks with VSL < 5 juim s~^ or 
VCL < 25 [xm s~^ (cut-offs imposed to avoid non-progressively 
motile sperm cells or groups). We then used the first 50 tracks 
(including both single sperm cells and sperm groups) recorded 
from each donor male in subsequent analyses for all but two 
males: in the case of one male of each species, fewer than 50 
tracks met our criteria (P. maniculatus male, n = 30 tracks; 
P. polionotus male, n = 27 tracks). Sperm group size was then 
subsequently counted for each track and verified on at least 
five different frames per track. 

We used two-factor (group size and donor male), two-tail 
ANOVAs to assess the effect of each factor on sperm average vel- 
ocity (VSL), speed (VCL) and linearity (VSL /VCL) within each 
species. After identifying the sperm aggregate size that achieved 
the greatest average velocity (n = 7 cells), we then compared the 



average velocity of seven-celled aggregates (the null) with the 
average velocity of all other sizes for each species using a one- 
sample two-tailed Mest. Next, we split the P. maniculatus and 
P. polionotus data into two groups and used a linear regression 
(with donor male as a covariate) to test the significant relation- 
ship between group size and average velocity at or below the 
optimum (n < 7 cells), and above the optimum (n > 8 cells). To 
identify how sperm aggregate size varies between species, we 
first averaged group size over each donor male, then used a 
two-sample two-tailed f-test to compare means, and an f-test 
to compare variances, of P. maniculatus and P. polionotus sperm 
aggregates. Finally, we used a two-way ANOVA (species and 
donor male) to compare difference between average linearity 
achieved by P. maniculatus and P. polionotus males. All statistical 
analyses were performed in R [38]. 

3. Results 

(a) Mathematical model 

We integrated equation (2.1) numerically for a wide range of 
parameter values. Our model sample consists of 100 cells in a 
square domain of size L = 500 (in units of the particle minor 
semiaxis b) with periodic boundary. For all choices of par- 
ameters, aggregation always leads to a prominent increase 
in the average velocity, VSL (but not speed, VCL), for small 
aggregate size, whereas large aggregates suffer from both 
reduced velocity and speed (figure 3d-i). The origin of this 
behaviour can be explained by noting that sperm can associ- 
ate with each other via soft adhesive bonds, modelled here as 
finitely extensible springs (see §2a). Once they are linked, they 
form aggregates whose structure is predominantly dictated 
by the geometry and the spatial distribution of the adhesive 
patches. The structure of the aggregates affects how the vel- 
ocity of the individual sperm determines the final velocity of 
the aggregates. Thus, radially symmetric aggregates consisting 
of many sperm (e.g. figure 3c) are likely to be non-motile, 
because the velocities of the individual cells effectively cancel 
each other. SmaUer aggregates, on the other hand, are asym- 
metric and maintain the typical head /tail directionality of 
individual sperm (e.g. figure 3b). More importantly, their 
close packed structure reduces the random fluctuations in the 
swimming direction of the individual cells, resulting in a per- 
sistent linearity of the trajectory and therefore a higher 
average velocity (VSL). The combined effect of these two 
competing mechanisms leads to an optimal aggregate size. 

The precise value of the optimal aggregate size, as well as 
the sharpness of the velocity peak, depends on the detailed 
geometry of the head /mid-piece complex and the adhesion 
properties of the sperm heads. To investigate how cell geome- 
try affects the swimming performance of an aggregate, we 
simulated self-propelled particles of various aspect ratios, 
the ratio between the length of the major and minor semiaxes 
of the elliptical particle. Increasing the slenderness of the par- 
ticles moves the velocity optimum towards larger aggregate 
sizes and simultaneously reduces the slope of the speed 
curve (figure 3d,g). This is because slender elliptical particles 
can pack more tightly than circles in two dimensions, so 
that it requires a larger number of particles to reach a sym- 
metric conformation. Increasing the adhesion range fa (thus 
the stall force that a single adhesive bond can withstand) 
also has the effect of moving the optimum towards larger 
aggregates (figure 3e), while leaving the speed essentially 
unaltered (figure 3h). Increasing the sperm propulsion 



Table 1. Two-factor two-tailed ANOVAs on sperm performance data. 



measure 


species 


factor 


F 


d.f. 


p-value 


average velocity (VSL) 


P. maniculatus 


group size 




1 


y.4l X lU 






donor male 


Id./ 


Q 
0 


< 1 X lU 




P. polionotus 


group size 


4z. 1 


1 


z.jz X lU 






donor male 


/.O 


Q 
0 


i.yu X lu 


speed iVLLj 


P. maniculatus 


group size 


z. 1 


1 


U. Id 






donor male 


r 7 
J./ 


Q 
0 


O.OZ X lU 




r. poiionoius 


group size 


2 ^. 
J.J 


1 
1 


U.UO 






donor male 


4.9 


8 


8.41 X 10"^ 


linearity (VSL/VCL) 


P maniculatus 


group size 


89.4 


1 


<1 X 10"^^ 






donor male 


17.4 


8 


<1 X 10"^^ 




P polionotus 


group size 


131.4 


1 


<1 X 10"^^ 






donor male 


16.5 


8 


<1 X 10"^^ 



velocity Vq (figure 3/) affects the position of the optimum 
only slightly, but produces a substantial improvement in 
the dynamic performance of aggregate. 

Finally, we note that in our two-dimensional self-propelled 
particles model, the aggregate size at which the speed starts to 
drop has a straightforward geometric interpretation related 
to the kissing number of the particles, defined as the number 
of particles that can touch a given central particle without over- 
lap. If the ellipses are not excessively slender, then this number 
equals six (the same as for circles), thus aggregates formed by 
six or more ellipses tend to be highly symmetric and undergo a 
severe drop in speed (figure 3g-i). Peromyscus sperm cells have 
a flat head-shape roughly similar to a very oblate ellipsoid 
(figure lb,c). For this type of shape, one might expect a kiss- 
ing number between six and twelve, the latter being the 
kissing number for spheres in three dimensions. 

(b) Experiments 

In both P. maniculatus and P. polionotus, motile sperm groups 
varied in size, ranging from 1 to 35 cells per group. We found 
that, overall, group size significantly influences average vel- 
ocity (VSL) in both species, even after the variation between 
donor males is accounted for (table 1). However, there is no 
significant relationship between the number of sperm in a 
group and speed (VCL), yet, similar to the result for average 
velocity, we found a significant effect of donor male in both 
species for speed (table 1). Finally, when we measured the 
linearity (VSL/VCL) of all sperm groups, we found that 
group size significantly affects linearity, with donor male as 
a covariate, in both species (table 1). 

The greatest average velocity was achieved by groups of 
seven sperm cells (figure 4a) and aggregates both smaller 
(n < 6 cells) and larger (n > 8 cells) than this number were 
slower in both species (P. maniculatus f = 4.2, d.f. = 8, p = 
0.003; P. polionotus t = 10.4, d.f. = S, p = 0.0001). Moreover, 
we found a significant positive association between sperm 
aggregate size and average velocity in both species as 
group size increased from one to seven sperm cells (P. mani- 
culatus = 0.48, p = 2.0x 10"^^; P. polionotus = 0.42, p = 
2.0 X 10 ~^^), yet a significant decrease as groups grew larger 
than seven cells (P. maniculatus R^ = 0.39, p = 1.63 x 10~^°; 
P polionotus R^ = 0.14, p = 5.20 x 10"^). 




p. maniculatus P. polionotus P. maniculatus P. polionotus 



Figure 4. Size and performance of Peromyscus sperm aggregates. 
{a) Mean + standard error (s.e.) of velocity (VSL) of sperm aggregates by 
group size over all donor males from each species, (b) Mean + s.e. group 
size of aggregated sperm in each species; sperm from each donor is rep- 
resented as a separate point with error bars, (c) Mean + s.e. linearity 
(VSL/VCL) of aggregated sperm over all males from each species; note 
truncated y-axis. (Online version in colour.) 

When we averaged sperm performance for each male 
in both species, we found that mean aggregate size did 
not differ significantly between species (mean + s.e. = cells 
per group: P. maniculatus = 6.0 + 0.2, range = 2-26 cells, 
P. polionotus = 6.5 ± 0.72, range = 2-31 cells, n = 50 groups 
per male, n = 9 males; p = 0.51, d.f. = 8); however, the aver- 
age group size in P. polionotus is significantly more variable 
than in P. maniculatus (figure 4h; P-test: p = 0.044). Moreover, 
the average linearity (VSL/VCL) achieved by sperm of 



p. maniculatus males is significantly greater than the average 
linearity of P. polionotus sperm (figure 4c; F = 47 A5, d.f. = 1, 
p=l.n X 10"^^). 



4. Discussion 

Our combined theoretical and experimental approach 
allowed us to build a mathematical model based on biologi- 
cal observations with testable predictions. Our empirical 
results are consistent with the salient predictions. These 
include: (i) when sperm conjugate in a head-to-head for- 
mation, such as in Peromyscus, and when aggregate size 
exceeds the optimum, cells will exert opposing forces upon 
one another and thereby reduce the velocity of the entire 
group, and (ii) the optimal size is dictated largely by the geo- 
metry of the sperm heads, and therefore species with similar 
sperm heads, such as P. maniculatus and P. polionotus [25], 
will achieve the same optima. We also found that, overall, 
group size significantly influences average velocity (VSL) in 
both species, and the greatest average velocity is achieved 
by groups of seven sperm cells — aggregates both smaller 
(n < 6 cells) and larger (n > 8 cells) than this number are pro- 
gressively slower in both species. Taken together, our results 
suggest that the shared aggregate geometry of P. maniculatus 
and P. polionotus sperm (probably as a result of the similarly 
shaped (figure lh,c) and sized [25] sperm heads) results in a 
similar relationship between sperm group size and 
performance, and thus similar optima, in these species. 

A third prediction of the model is that sperm aggregates 
achieve greater average velocity (VSL) because they move 
in a more linear path, rather than an increase in speed 
(VCL). Indeed, we found no significant relationship between 
the number of sperm in a group and speed in either 
P. maniculatus or P. polionotus. However, like velocity, we 
found a significant effect of donor male in both species on 
speed; this variation among males and between species is 
consistent with earlier findings of intermale differences in 
speed in these species (Fisher et al. 2014, unpublished data). 
By contrast, when we measured the linearity (VSL /VCL) of 
all sperm groups, we found a significant effect of group 
size on linearity with donor male as a covariate in both 
species. These results indicate that the benefit of sperm aggre- 
gation is, indeed, conferred via a more direct path of travel, 
rather than a change in speed, as predicted by our model. 

Our experimental results are consistent with the predic- 
tions of the model that shared aggregate geometry of 
P. maniculatus and P. polionotus will yield similar relationships 
between sperm group size and performance in both species. 
In nature, however, the ideal strategy is not always the 
most prevalent one owing to associated costs, selection on 
pleiotropic traits and/or genetic drift. While P. maniculatus 
sperm have evolved under a selective regime with intense 
competition [21,22], evidence suggests that monogamy in 
P. polionotus [23,24] is derived [39] and, therefore, sexual 
selection is likely relaxed in P. polionotus. When we measured 
the average aggregate size in each male across the two 
species, we found that while the average group size does 
not differ significantly between species, the average group 
size in P. polionotus is significantly more variable than 
P. maniculatus. In other words, the mean group size does 
not differ between the species, which are both within one 
cell of the observed (empirical) and predicted (theoretical) 



optimum, but the distribution around the mean is signifi- 
cantly larger in P. polionotus, and thus more aggregates 
are further away from the optimum, compared with 
P. maniculatus. These results suggest that sexual selection, 
and male -male competition specifically, may be imposing 
stabilizing selection on aggregate size in P. maniculatus 
sperm, resulting in more groups at or near optimal size; by 
contrast, the monogamous mating system of P. polionotus rep- 
resents a relaxation of male competition and is consistent 
with greater variation in sperm group size. 

Given that P. maniculatus sperm are more likely to 
form aggregates at or near the optimal size compared with 
P. polionotus, our model also predicts an overall increased 
linearity in P. maniculatus sperm in the total sample. 
Indeed, P. maniculatus sperm move in a more direct trajectory 
(VSL /VCL) than P. polionotus sperm. The results from this 
study reveal that selection may, therefore, act on sperm 
swimming performance via aggregation behaviour. 



5. Conclusion 

Our detailed observations of sperm shape, aggregate geometry 
and their dynamical performance suggest an optimal sperm 
aggregate size that leads to a maximum linear velocity of a 
group. Our minimal mathematical model — that accounts for 
the geometry of the sperm and the mechanics of their adhesive 
interactions, when combined with the d3mamics and fluc- 
tuations of translational and rotational movement — captures 
the non-monotonic dependence of aggregate velocity on the 
number of sperm in a group. The underlying mechanism is 
simple: in small groups, sperm adhesion reduces the size of 
rotational fluctuations by effectively cancelling them, whereas 
in large aggregates, this effect is dominated eventually by redu- 
cing the mean translational velocity owing to the isotropic 
geometry of a cluster. Thus, relatively few mechanical par- 
ameters can explain the dynamics of a seemingly complex 
biological process. 

Our empirical data test these model predictions and show 
that sperm achieve greater velocity surprisingly not by 
increasing speed, but rather by travelling in a more direct 
path than solitary cells. This collective behaviour arises 
from direct physical interaction among cells, which deter- 
mines the optimal aggregate size. The number of cells 
involved in an aggregate, therefore, greatly contributes to 
sperm performance and the reproductive success of a male, 
thereby offering another dimension of sperm biology on 
which selection can act. Moreover, by comparing sperm 
dynamics in two species that have evolved under disparate 
competitive regimes, we are able to implicate a role for 
sexual selection in the evolution of complex behaviour of 
these seemingly simple cells. 

Thus, a deep understanding of sperm behaviour requires 
us to combine our knowledge of geometric and physical 
constraints with reproductive biology; indeed, these 
dynamics are clearly driven by a combination of morphology, 
kinematics and the selective environment. While selection 
ultimately acts on organismal fitness, our picture allows us 
to link this to the dynamics of movement and the adhesive 
interactions among sperm. Indeed, sperm cooperation and 
competition are a remarkable arena not only to study evol- 
utionary processes, but also to test quantitative models for 
how they may play out in nature. 
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Appendix A 

(a) Single sperm dynamics 

To understand why aggregation improves the mobility of 
sperm, it is useful to consider the dynamics of a single cell 
resulting from equation (2.1) in the main text: 

^ ' (Al) 



dt 



-■ Voti and 



dt 



where ^ is a zero-mean delta-correlated Gaussian random 
variable: U{t)^{t')) = 2D8{t - t'), with D a rotational diffusion 
constant. Equation (A 1) describes the motion of a cell moving 
at constant speed Vq, but whose direction of motion is affected 
by random rotations. Translational noise can be neglected 
due to the fact that the typical Peclect number of a swimming 
sperm of size L is Pe = VqL/Di = 10^ — 10^ [40], where Dt is 
the translational diffusion constant of the sperm; thus, diffu- 
sion is negligible compared with drift. The beating of the tail 
that controls the orientation of the cell, on the other hand, is 
subjected to fluctuations owing to the noise in the activity of 
the motors regulating the flagellar beating. Equation (A 1) can 
be easily solved using the standard machinery of Brownian 
motion [41,42]. Integrating the 0 equation yields 



AO{t) = e{t) - 0(0) : 



from which 



{Ae(t)} = 0, {A0(t)A0(t')} = 2D mm{t, t'), 



(A 2) 



(A3) 



where min(f, t') represents the minimum between the times 
t and t' . Now, as A 0 is a linear combination of random 
variable with zero mean and variance {Ae\t)) = IDt, from 
the central limit theorem, it follows that it is Gaussianly 
distributed. Namely 

1 



3-(A02(O)/(4DO 



(A4) 



where we called P[/^0{t)] the probability density function 
associated with A 0(f). This allows us to calculate the averages 
of exponential and trigonometric functions 



(e 



±iA0(O 



3-(l/2)(A02(O) ^ 



from which 

( cos 0{t)) = e"^^ cos 00, ( sin 0(t)} = e"^^ sin ( 



(A5) 



(A 6) 



Averaging both sides of the position equation (A 1) allows us 
to calculate the average velocity of a sperm cell: 



{V(t)} : 



(AT) 



where n(0) = (cos0o/ sin0o) is the direction of the cell axis at f = 0. 
Equation (A 7) represents an important property of our model 
sperm cell: at short times, individual sperm move balHsticaUy 
in the direction of their axis, but, owing to the noisy flagellar 
beating, their velocity exponentially loses directional 



correlation and in a time scale of order = 1 / D the cell has com- 
pletely lost track of its initial direction. Aggregation allows 
groups of cells to reduce the fluctuations in the direction of 
motion, thus increasing the correlation time (see later). 

The mean-square displacement of a single cell can be 
calculated as 



{\r(t)-riO)\'} = vl 
and yields 



{cos[0{h) - 0{t2)]} = e 



dtidt2{cos[e{ti) - 0(^2)]), 



(A 8) 



-{l/2){[^^{h)-^^{t2)f) 

-D[ti+t2-2 min(ii,i2)] - 



-D\h-t2\ 



Combining this with equation (A 8) gives 

^-Dt\ 



{\r(t)-r{Or}- 



D 



(A 9) 



(A 10) 



Thus, for Tc, the motion is completely Brownian with 
{\r{t)-r{0)\^)'^lvQt/D. The average speed (VCL) and vel- 
ocity (VSL; figure 2) of a single cell are then given as follows 



VCL = i7o and VSL: 



{\r{t)-r{m 



(All) 



Because the mean-square displacement does not depend lin- 
early on time, VSL depends in general on the time range used 
to calculate it and vanishes in the limit f ^ oo, indicating that 
a cell has been diffusing long enough to visit its initial pos- 
ition more than once. For the dynamics of a single sperm 
cell, we can replace t in equation (A 9) with the lifetime ti 
of an isolated cell, defined as the time it takes for the cell to 
meet other cells and form and aggregate. If the sample con- 
sists of a density p of uniformly distributed cells, then this 
is approximately given by ti ^ \/{Voy/p). Figure 5 shows 
the average velocity and the lifetime of isolated cells for var- 
ious values of Vq and cell number obtained from numerical 
integration of equation (A 1). 

In the experiments, the time range ^exp used to calculate VSL 
is fixed and is based on the frame rate of the images (30 frames 
per second). This can be incorporated in the definition of D, by 
introducing a dimensionless rotational diffusion number 
V = Df exp • The long-time average velocity is then characterized 
by two quantities: VCL = Vq and VSL/ VCL ^ 1/VV. Further- 
more, the above-presented discussion is not only restricted to 
the case of an isolated cell, but can also be extended to aggre- 
gates. In the case of an aggregate, VCL is the average velocity 
of the aggregate, whereas VSL /VCL can be taken as a measure 
of the effective diffusion number Peff that accounts also for the 
reduction in the orientational noise resulting from aggregation. 



(b) Scaling 

The scaling behaviour of the VSL for increasing aggregate 
size can be predicted within a mean-field framework starting 
from a simple geometric argument. As explained in the pre- 
vious section: VSL ^ ^ 1/AO. Thus, if the angular 
span A0 is reduced by a factor x < 1, then VSL is increased 
by a factor 1/x. Now, in the interior of an aggregate, A0 
decreases owing to the crowded environment. To account 
for this effect, we use the following ansatz: 



(A 12) 




^0 ^0 

Figure 5. Average velocity (VSL) and lifetime (i.e. aggregation time) of a single cell. Both VSL and Vq are expressed in units of q/tc, with a the length of the cell 
minor semiaxis and = MD the correlation time. The circles are obtained by numerical integration of equation (A 1), whereas the solid lines are given by equation 
(A 11) by setting t = A/(vq^), with p the particle density and A a fitting parameter. 



10 



> 10~ 



10- 











^0 






A 1.5 






□ 1.0 






o 0.5 







5 



10 



20 



Figure 6. Log -log plot of the average velocity (VSL) versus aggregate size for 
various values of propulsion velocity Vq. The numerical data are the same as in 
figure 3f, while the dashed line shows the predicted scaling law VSL 
two dimensions. (Online version in colour.) 



where (2^ is the solid angle in ii-dimensions, the ii-dimensional 
kissing number and n the average number of neighbours 
in the aggregate. Equation (A 12) implies that each new 
neighbour takes an equal amount of angular space until, for 
n = kd, there is no space left and the cells have reached a 
jammed configuration. 

The number of neighbours n depends in general on the 
size N of the aggregate and in particular on the number of 



cells in the interior of the aggregate compared with those dis- 
tributed along the boundary. Thus, calling and the 
average number of neighbours in the interior and on the 
boundary, we have 



ZjNj + ZbNb 
N 



(A 13) 



where Ni and Nb represent, respectively, the number of cells 
in the interior and on the boundary and N = Ni + Nb. Next, 
assuming = k^ and taking into account that Nb ~ N^^*^, 
we obtain 



-cfc-Zb)N-(^-i)/^ 



(A 14) 



with c a constant. Now, when the aggregate consists of a single 
cell, n = 0. This allows us to calculate the constant c in the 
form c = kd/{kd- Zb), from which n = k^il - N~^^~'^'>^^) and 
equation (A 12) becomes 



A^-/2rfN-^^-i^/^. 
Consequently 

VSL(N) = VSL(l) N^^-i)/^ ^ 



JnV^ 

\n2/3 



(A 15) 



(A 16) 



For our simple two-dimensional model, this prediction is 
consistent with the numerical data shown in figure 6. 
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